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ABSTRACT 

We construct the classical r-matrix structure for the Lax formulation of BCat Ruijsenaars- 



Schneider systems proposed in [|18|. The r-matrix structure takes a quadratic form 
similar to the A^r Ruijsenaars-Schneider Poisson bracket behavior, although the dy- 
namical dependence is more complicated. Commuting Hamiltonians stemming from 
the BCat Ruijsenaars-Schneider Lax matrix are shown to be linear combinations of 
particular Koornwinder-van Diejen "external fields" Ruijsenaars-Schneider models, 
for specific values of the exponential one-body couplings. Uniqueness of such com- 
muting Hamiltonians is established once the first of them and the general analytic 
structure are given. 
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1 Introduction 



The relativistic extensions of Calogero-Moser A^-body integrable systems, originally 
introduced by Ruijsenaars and Schneider have been the subject of numerous 
investigations in these last years. Their exact connection to field-theoretical integrable 
systems, initially described in was clarified in P]; their dynamical classical r- 
matrix structure, first tackled in 0, Q was finally established in and characterized 
as a quadratic structure, a la Sklyanin |P, 0, stemming from the dynamical linear 
r-matrix structure of Calogero-Moser systems 0. This in turn is connected to the 
realization of such systems as "Hamiltonian reduction", in a more general sense 
see |T^, of dynamical systems living on a Heisenberg double |T^, |1T|, where the 
quadratic r-matrix structure is natural ||10[- More complete descriptions may be 



found in |]T3| . 

More recently several extensions of the Ruijsenaars-Schneider systems were con 
sidered. First of all, one-body-potential (or "external field extensions") were con 
structed in 
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The integrability proof, and construction of the quantum Hamiltoni- 
ans themselves, rested upon polynomial-algebraic arguments pioneered by Koorn- 
winder and indicating a connection to BCTv-type algebras. Both Lax formulation 
and classical r-matrix construction were lacking. Quantum elliptic van Diejen-type 

T7[ | using corner-transfer-matrix methods. 



Hamiltonians were then constructed in 
and incidentally pointing again to a BCat 
ticular Koornwinder-van Diejen potentials, 
in 



structure underlying at least some par- 
The corner-transfer-matrix method used 



17] however does not exhibit a clear-cut way of defining a classical limit with one 



single Lax matrix, using as it does two types of Lax operators. Finally an explicit 
construction of classical BCat and Cn Ruijsenaars-Schneider models was presented 
using a Z2-folding of the original A2Ar or A2iv+i algebra intrinsically present 



m 
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in original Ruijsenaars-Schneider models. The Lax matrix and classical commuting 
Hamiltonians were then constructed, again as Z2-foldings of the A2Ar (A2Ar+i) Lax 
matrix and commuting Hamiltonians. The elliptic BCtv and systems were then 
defined in 



IS 
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and the D^r case was considered in 
This series of results still begs several questions and we wish to address here two 
points which remained unclarified. 

First of all we construct the classical r-matrix for hyperbolic BCat models (we 
expect that the Cjv and D^r case, and the elliptic potentials, may be treated by 
similar techniques although the elliptic case may endow more complicated algebraic 
manipulations). This problem may seem academic, since in any case commutation of 

However it actually sheds 



18 



the Hamiltonians follows from the construction itself 
light on the delicate question of interplay between the folding procedure A2Ar — > BCn 
and the initial "quadratization" of Poisson structure entailed by the change of base 
symplectic manifold for the "Hamiltonian" reduction from T* g (cotangent bundle Lie 
group) to D qh (Heisenberg double of Lie group) [[ Indeed it will eventually turn out 

"'^ Similar difficulties occur when considering the quantum deformation of BC- or D-type algebras, 

compared with deformation of A^r algebras. This fact was pointed out to us by D. Arnaudon. 
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that the classical r- matrix for hyperbolic Ruijsenaars-Schneider BCa? models may 
be recast under a quadratic form, similar to the A/^r case, but contrary to what 
occurred in the Calogero-Moser models |21| the structure (after folding) exhibits now 
a dependance in both sets of conjugated dynamical variables, and the quadratic r- 
matrix is thus not directly related to the linear dynamical r- matrix structure for BCiv 
Calogero-Moser models. 

The second problem which we consider here deals with the connection between 
the initial Koornwinder-van Diejen Hamiltonians and the "canonical" Poisson com- 
muting Hamiltonians generated by the traces of powers of the Lax matrix for BCn 
systems. It will be shown that the Koornwinder-van Diejen Hamiltonians are in fact 
combinations of the "canonical ones" , and this property is actually true for any set of 
Poisson-commuting Hamiltonians with the same functional structure (to be explicited 
hereafter) . 



2 The classical r-matrix structure 



2.1 The BCiv Ruijsenaars-Schneider models and notations 



The canonical variables are a set of rapidities {6i, i = 1 ■ ■ ■ N} and conjugate positions 
Qi such that {9i,qj} = 6ij. The Hamiltonian reads: 



H 



N 

E 

j = 1 

e = ±1 



-£(31 



(1) 



where 



1 1/2 



N 



filj) /(2 Qj) Yl f^^i ~ ^'^) f^^i + 



fc = 1, 



N 



and W = JJ f{qk). 



k=l 



Function / may take different forms, namely: 



/(g) 
fiQ) 
fiQ) 



(rational) 



sinh 7 
sinh^ Y 



(hyperbolic) 



sin^7 

1 7T—r (trigonometric). 

sin ^ 



The most general elliptic case where: 

/(g) = (A + vV^q)) , V = Weierstrass function 
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will not yet be considered here. 

Since the trigonometric and hyperbolic cases define the same model at least locally 
up to a redefinition of the parameters 0, and the rational case is obtained by an easy 
limit procedure from one of these models, we shall consider in the following only the 
hyperbolic version. 

Let us note that one can also write /(g) = v{q) v{—q), with : 



v{q) 



sinh(^ + 7) 
sinhy 



or even as a rational fonction of an exponential variable: 

This rational formulation will be useful to establish some functional identities a la 
Liouville 0. 



2.2 The BCat Lax operator 



As shown in [18|, the Lax formulation of BCtv Ruijsenaars-Schneider system may be 
obtained as a folding of the A2„ case. The reduction works as follow: labelling the 
2 + 1 rapidities {9i, i = —N ■ ■ ■ N} and conjugate positions {g,, i = —N ■ ■ ■ N}, one 

sets 9i = Si 6\i\ and qi = Si q\i\ with: 

+1 for 1 < i < A^ 
for i = 
-1 for -l> i> -N 



The Lax matrix for the A^^„ cases is: 



^2n 

N 



i,j=-N 

N 



Cij{qi,...,qN,qo,q-i,...,q-N,Oj) = c{qi-qj)e^^' JJ f^^iqj-Qk) (2) 



k = -AT, 



where {etj} is the standard basis for (2 A^ + 1) x (2 A^ + 1) matrices; / was given in 
the previous subsection and 

e(,) = = (1 - A) 



sinh(^ + 7) ' ' z-X 



^ The global structure of trigonometric vs hyperbolic models is however quite different, owing to 
topological properties, as can be seen for instance in p3| 
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The Lax matrix for the BCat Ruijsenaars-Schneider systems then reads: 

N 

L= ^ LijCij with Lij = Cij{qi, ...,qN,0,-qi, ...,-qN,ej9\j\) (3) 

i,j=-N 

It can be rewritten: L^j = c{qi — qj) e"^^-''^'^' fj, extending the definition of fj given 
in (P to J e {-A^...A^}. 

Note that with this extension of fj one has fo =U and fj = f-j. 

It has been shown that the Lax operator (|) satisfies the quadratic fundamental 
Poisson bracket 0: 

{Cf C} = C0Cai-a2C0C + C2SiCi-CiS2 A, (4) 

where Ci = C ® 1, C2 = 1 ® C and the quadratic structure coefficients read: 

ai = a + w, Si = s — w, 

a2 = a + s — s'" — w, S2 = s'^ + w . 

For any matrix M, the matrix M'^ is defined by: 

N N 

if M = ^ Mijki Cij ® Cki then M'" = ^ Mijki eu ® Cij. 

ijkl=-N ijkl=-N 

Matrices a, s, w take the form: 

AT ^ 
a = a ^ coth-(gfc - gj) Cjfc (g) Cfcj, 

j, k = -N 



N 

-a 



N 

V 

w = a 



coth-{qk - qj) Ckk (S) Cjj (5) 

J, fc = -JV 

V 

where a = /3 — . 

It must be recalled here that the most general structure of Poisson bracket for a 
Lax operator of a Liouville-integrable system is a linear one [^: 

{L?L} = [r,Ll]-[r^L2]. (6) 

The quadratic form corresponds to the case where the r-matrix itself assumes a 
linear dependency in L of type: 

r = bL2 + L2C, (7) 
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with b and c arbitrary matrices determining the quadratic coefficients ai, a2, si, S2.: 

ai = — c, a2 = b'^ — b, si = c + and $2 = s\. 

In the next subsection, we will show that the BCat Ruijsenaars-Schneider Lax op- 
erator (I) also satisfies a quadratic fundamental Poisson bracket albeit with a 
fundamental difference with respect to regarding the dependence on the dynam- 
ical variables. We will give explicitely the generalizations of the matrices ai, 02, si 
and S2, hereafter denoted "quadratic r-matrices" for obvious semantic reasons. 



2.3 Computation of the classical r-matrix 

Let us calculate the Poisson brackets of the Lax matrix @: 



rr T ^ pt t f 9\nLij d\nLki\ 
V dqyi J 



and express it in terms of the Lax matrix (|^) 



d In Lij / (9 In Cij d In Cij 



dq\i\ V dqi dq^i 
We thus get: 

rr r 1 _ r/- r \ ^ ar t f^^^^ki dlnCij\ 
i-i^ij, J^klf — i^ij,^klf -r P J^ijJ^kl y j ■ 

The Poisson bracket of the first term on r.h.s. keeps the same form (^) where one 
should fold the dynamical variables {6i = Si 6\i\ and qi = Si q\i\). We thus only need to 
concentrate on the remaining term, introducing the four-index object: 

^ 2 / d In Cki _ d In dj 

'-^ijkl — I o o 

u \ dq^j dq^i 
Straightforward calculations yield: 

Uijki = Sj-iUj + {Si -I — Uj — {6j-k — Sj-i) tki 

where 

tij = -- (Inc)'(gi - qj) = coth(^(gi - qj) + 7) = ^'^^^^ 



2 ^ 



= — 

V 



k = -JV, 
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N 

2 coth-(gj - qk) + coth(-(gfc - qj) + 7) - coth(-(gj - qk) + 7) 



fc = -AT, 



^ 7 T + 7 ~ T VV ~ 1^ ^ + - 

— Zyfc Zk — A Zj Zj — AZk 

k = -N, k = -JV, 



fc # i fc / i 



with ajfc = coth— (gj — qk) = ^ for j 7^ fc. 

Z Zj Zk 

Note the following properties of these objects on the folded space: 

In the expression of Uijki, the terms in 6j^-i are not on the same footing as the 
others since they are separately antisymmetric under the exchange of the two spaces 
(operation tt) whereas the remaining terms verify this property only altogether. 
We thus first take care of the 6j-i terms, introducing the matrix p: 

^ 1 
p = a ^ Lki{-ui- tki) e^i^i (g) eki, 

k,l=-N 

realizing them linear r-matrix form (0): 

{[p, Li] - L2])ijki = a LijLki (^Sj^^i {uj - Uj + tki) + ui - tki) - Sj^-k Uj - Uj)^ 

= a LijLki {Uijki — Uijki), with Uijki = Si^^i {tij + tki — - ui) — Sj^^k {tij +tki — - Uj). 

We may furthermore bring it back to our seeked-for general quadratic form by setting 
p = T L2 (i.e. taking h = r and c = in (0)), since the matrix L is invertible: 

T = pL2^=a ^ Lk-iLz]i{-U-i-tk-i)eii®eki. (8) 

i,k,l=-N 

One should immediately note, from the explicit form of L (^), that this matrix r 
actually does not depend on the rapidities ^j's. We are therefore still in the "canoni- 
cal" quadratic structure (a la Suris) of type (^ with dynamical quadratic r-matrices 
depending only on one set of canonical variables (the ^j's). 

We will now show that there exists a matrix cr, such that cr'^ = a and: 

LijLki Uijki = [L2 cr Li — Li a'^ L^l^jki ■ 

This corresponds to setting c = a and 6 = in (0) and thus actually formally 
completes the quadratic r-matrix structure (H). 
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In order to ensure self-consistency of the dependence in the indices, we assume 
the following tensorial structure: 



TV 

(j= ^ amnemn® e-n-m, Satisfying cr" = cr ie: a^n-m- ^mn, 

m,n=—N 

yielding: 

N N 

[L2 a Li — Li a'" -^2)1^^; = ^i-l CFin Lk-n Lnj — ^j-k CTfcn -^i-n Lnl- 

n=-N n=-N 

The set of equations to be solved then reads: 
or, equivalently: 



^ T T ■ 1 



V] CTm " - Oi {tij + + - lij) = (5j Sjfe, with Sjfc = Ski- 

J-'ijJ-'k-i ^ 



n=-N 

Direct calculations yield: 



T , T pl^^" f/ f . — f , ■ — 1 

' (1 - k-k) -r — 7 — + ^j,-k -i 



LijLk—i g/^^i Y ti—k tj^j 1 J 

and equations become: 

^ 1 1 

^ ^ ^in i^nj tn—k) — ^ {tij ^i— fe 2 ^^'^ 2 ^i—k) 

n=-N 



t ■ — 1 1 



n=-N 

Let us recall here that the only additional contraints on matrices a and s read: 

_ ff ~ 

" 11 (9) 

Sij — Sji 

According to the previous equations, the matrix s is determinated given the matrix 
a, and a is obtained, up to a one-dimensional degree of freedom Vi, by: 

^ 1 

^nj = a {tij +2'^' + ^0 ' (10) 

n=-N 
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since t is invertible. 

Remains to verify that one can find f j's, such that the compatibihty relations @ 
be satisfied, namely: 

N N 

^ ^ tim fm ^mj (tmj tim ~l~ ^m) ^ ^ fm (^— m tmj ^jm) (-^-^) 

m=—N m=—N 
N N 

(^^ + 1) E (^-.™ - 1) r ' w = (^-. + 1) E (^L-i)r^) 

n,m=— Af 



(12) 



n,m=—N 



N 



Equation directly yields Vi = —1 + r] ^ (t^?^ — 1) (t with 77 an arbitrary 

n,m=—N 

constant. 

In order to solve (0) we shall first compute its left hand side. 
We do so, by twofold evaluation of the following contour integral in the complex 
plane around infinity: 



I 4 n 



2m 



Co, 



dz Zi + Xz z + Xzj 
z Zi — X z z — ^ 



N 

n 

k=-N 



Z — XZk Zk — Xz 



z - Zk X{zk- z) 



Contour Coo is oriented counterclockwise and loops around infinity. 
Residue at infinity gives lij = — 1. 

Whereas lij also equals the sum of residues at poles of the meromorph integrand 
in the whole complex plane, that is: a single pole at z = with residue —1 and a set 
of double poles aX z = Zm- We thus obtain: 



N 

E 

m=-N 



Zi -\- X Z Z -\- X Zj {^Z — A 2^m) (^^m — A 2^) 
Zi — X Z Z — X Zj 



-Xz 



N 

n 

fc = -JV, 



z — A Zk Zk — X z 



z - Zk X{zk- z) 



, Z Z^i 



Noticing that: 



2 ^ 2 d ^ 

E (ln/)'(9fc-gm) = — o — In W f{qk-qv 



V 



k = -N, 
k ^ m 



2 Zn 



_d_ 

dZrr 



N 



m n 



Zm ~ X Zk Zk — Xz. 



k = -JV, 

k ^ m 



k = -JV, 
k ^ m 



Zk X (yZk Zmt) 



~2 Zri 



N 



i„ n 



z — Xzk Zk — Xz 



z- Zk Xizk- z) 



k = -JV, 
k ^ m 



[z = z„ 
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and also: 



, .Zi + Xz z + Xzj^ 
ln( 7 

Zi — A Z Z — X Zj 



we obtain: 



TV 



2 t ■ 



m=-N 



or equivalently: 



N 



N 



^ ^ iim fm ^mji^mj ^im ~l~ ^m) ^ ^ /^m ('^im i 



m=-N 



m=-N 



This derivation of Liouville-type functional identities stems from similar derivations 
to be found in IPI. 



Equation ([TTD now reads: 



TV 



TV 



N 



^ ^ ^im fm (^m ~l" 1) ^ ^ ('^— m ~l" 1) fm ^mj ^ ^ t—jm fm (^m ~l" 1); 



m=-N 



m=-N 



m=-N 
TV 



and can be solved straightforwardly: Vi = —1 + ^ {t ^)im-, with ^ 

J i 



any con- 



m=-N 



stant. 

We thereby prove the consistency of (§) and ([T0|) , since Vi = —1 is an obvious 
solution (with ?7 = ^ = 0). 

In addition, performing calculations of the same type of contour integrals, one 
gets: 

N 4A 1 ^ 1 

n,m=-N ^ ' m=-N « 



with D; 



AT 

n 

fe = -JV, 



Zi — X Zk 



The two forms which the fj's should satisfy, are actually identical: fj = — 1 + -^^7, 
C being an arbitrary constant. 

We fix this remaining gauge, setting ( = , in order to obtain the simplest 



1- A 



form for a from relation (fi^) - 



TV 



= " ^ Si - {I- Sij) aij) , where Si = j^-^ + ^ + 



m=-N 
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We finally give the expression of the matrix cr: 



^ A- [d+ „ 

cr = a ^ -J- {6ij Si - {1- 6ij) ttij) eij^e_j_i, with Ai = J e~^^\ (13) 

iJ=—N * V ~* 

The r-matrix structure is now completely defined by a quadratic Poisson bracket 
of type (H) where the quadratic r-matrices ai, 02, si and S2 are changed into: 

Oi — ai = ai, 

si ^ si = si + a + r'' 

S2 ^ S2 = Si = S2 + a"" + T = S2 + a + T 

a2 0.2 = ai + si - S2 = a2 + t"" - r, (14) 

and matrices a and r are respectively defined by equations (p!3D and (H). 



2.4 Comments 

One should notice that this quadratic r- matrix structure is now fully dynamical, 
depending both on the positions gj's and rapidities 6'j's. Moreover, its conjugating 
factor Ai, which bears this double dependance, is deeply connected to the structure 
of the matrix L under folding: 

r-i _ A' r . . 

^ij - -J- ^-^-3 

We have here an interesting first example of a "doubly dynamical" r-matrix depen- 
dence, stemming from the interplay between the folding procedure leading from A^v to 
BCat, and the quadratic structure of Ruijsenaars-Schneider-type models. This seems 
to open new perspectives, first of all on the classification of such doubly dynamical 
r-matrices. The only other example known to us at this time is the classical linear 
r-matrix structure for the Lax formulation of the A^r elliptic Calogero-Moser model 
in the absence of spectral parameter . Proof of its double dynamical dependance 



is given in although the explicit form is yet conjectural Curiously however, 
it occurs in relation with a,n An model with no folding and may therefore be of a 
different nature. 

Interpretation of doubly dynamical objects in the frame of quantum group theory 
is yet lacking. "Simply" dynamical r-matrices are known to be connected to the theory 
of Drinfel'd-twisted quantum groups, specifically of the type of Felder's Dynamical 
Quantum Groups (see for instance ^). Whether doubly dynamical objects 
have such connections is a new problem and we have no further comments to make 
on this point. A very recent result |^ may however give indications on how to twist 



quantum groups by non-abelian twisted cocycles (here the twist would occur "along" 
a Heisenberg algebra). 
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3 The canonical Hamiltonians 



3.1 Preliminaries 

We first describe the Poisson-commuting Hamiltonians generated by traces of powers 
of the BCtv Lax matrix (|^). They read for any integer / G {l..A^}: 

'Hi = tr{L')= Yl ^^(^)' 

J C {-JV..JV} 



with mj{L) the principal minor of L with lines and columns indexed by J^. 
Taking into account the form of L and properties of Cauchy matrices: 

mj{L) = e-^'^ n f"'^<i3-<ik). where Oj = Y,(^,. 

k!t J 



Because of the folding, we now rearrange these terms so as to sort them with respect 
to linearly independent exponentials of rapidities. We thus decompose J = eJ\JS, 
separating indices of J such that their opposite does not belong to J (set eJ = 
{£j\j\/j e J" A — j ^ J} and J = {\j\/j G eJ} C {l..A^}) and the complementary 
parts, symmetric under folding S: 

ni= Yl ^J^'HJi n /'^'(^.-^'^)' (15) 

J c {i..iv}, I J| <i j e sJ 

£j = ±1, j e J k ^ sj 

with 

^^'P= E n /'^'fe-?^)= E n <<is-<ik)- (16) 

S <Z Ak = K{j -K\j{0} seS S C Ak sgS 

s = -s, |S| = p fc s Ak\s s = -s, \s\ = p fee Ak\s 



We then recall the Koornwinder-van Diejen Hamiltonians in the classical case: 



Hi= E kYJ^ ^-S;^- (17) 



J C {l..iV}, I J| < ! 

Ej = ±1, j e J 



where, after some rearrangements. 



fc 6 ylx U -eJ 
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and 



TT ( ^\P ^ T[ ^f^^*) TT TT ( \ 

£-f C Ak i, i' S e/ i e el 

|/| = p i < i' fee -A/fVE-f 

The w are particular functions and may be interpreted as an interaction with some 
external field. 

Direct computation yields: 

v^-ej;j<= = TT , ^^^\ . ~r~J^\^ — V IT f^'ij ~ 

tk^j ^(2 qj) v{qj) v{-2 qj) v{-qj) J- J. 

Jfc'=>' j £ eJ 



Setting w{qj) = v{2qj) v{qj), which is an admissible choice according to []T3[, Hi (|T^ 
takes actually the same form as Hi (p!5D, up to the change of Uk,p in to Uk,p- Uk,p 
takes a simpler form, for this one body potential w: 

e/ C -4ft: i, i £ el i £ el 

\I\ = p i < i' k e AK\eI 

nevertheless it is not generally equal to Uk,p (the notation used in |]18| is in this respect 
misleading) . 

They are actually only equal for p = 0, where trivially: Uka = 1 = Uk^- For 
instance when p = 1, one gets: 

Uk,i = - ^ JJ v{qi - qk) and Uk,i = JJ v{qk). 

ie.AA'\{0} k€ Ak ke Ak 

k^i k ^0 

We compute a suitable contour integral on the same lines as in the previous section 
to obtain the Liouville-type functional identity: 

V- TT r ^ sinh7(2|ir| + l) 

z — / J. J. gmh 

teAK kGAK 
k i 

and thus show that: Uk,i = 1^k,i — ^"^^linh^^"*"^^ - 

It will now be shown that these two relations, for p = and p = 1, between the 
Uk,p's and Uk,is are actually sufficient to establish that the two sets of Hamiltonians 
define the same family of commuting dynamical fiows, namely one set of Hamiltonians 
is a triangular linear combination of the other set. 

A more general result will in fact be proved in the following subsection. 
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3.2 Uniqueness theorem 

Theorem 1 Let qi and 9i, i E N, be a set of conjugated variables such that {6i, qj} = 
6ij. Let I and K be arbitrary finite sets of indices included inN. Assume the existence 
of a set of complex functions uk,p depending upon the set of indices K and a natural 
integer p, and of another set of complex functions v^^jj depending upon the sets of 
indices J and I (J G I) and a \ J\-uple of signs e = {ej, j G J), such that: 

• uk,p and v^jj be independent of the rapidities 9iS. 

• Uk,o = 1, Vdtj = 1, and Vs{j}j ^ 0. 

• = { hj = e"^^""' v^jj, I E } be a family of Poisson- 

J C I, \ J\ < I 

£j = ±1, J e J 

commuting functions (O^j = ^j^j)- 

If there exists a second set of complex functions uk,p obeying the first two condi- 
tions; such that = { hj = ujoi_^j^ e~^^^-^ Vsjj, I G {l..|/|} } be a new 

J C I, \ J\ < I 

Ej = ±1, J e J 

family of Poisson-commuting functions; and Uk,i = Uk,i + Cid-fTl), then there exist 
coefficients Cr{m), (r, m) G N^, independent of all dynamical variables, connecting the 

two families of Hamiltonians as: 

I 

hj = Q_s(|/| — s) hi, with Wm G N, co(m) = 1. 

s=0 

Proof: 

The strategy of the proof rehes upon a recursive procedure on p, showing that: 

p 

uk,p = Cp_r{\K\ — r) UK,r, C N finite and such that \K\ > p. (19) 

r=0 

Let Iq be a strictly positive integer; the recursion hypothesis hereafter denoted r.h., 
states that (0) is valid for any p < Iq. 

The assumptions in the theorem immediately imply the validity of r.h. for Iq = 1, 
can be directly rewritten as: 

hi = hi + c,{\I\). (20) 
Let us assume r.h. up to Iq and establish it for Iq + 1. We have: 

J c I, 1 < \ J\ < lo + 1 
Ej = ±1, j e J 

13 



Since Iq + 1 — \J\ < Iq in the previous summation, we apply r.h. to get: 

J C I, 1 < \J\ < s 

= ±1, j e J 
lo+l lo+l 
= - - s)m7,, + ^Qo+l-.(|/| - S) /l^. (21) 

s=l s=l 

We now use the Poisson-commutation property of as: 

{hio+i,hi} = 0- 

Combining pO|), (|2T|) and the Poisson-commutation property of yields: 

lo+l 

= {^i,io+i - X] Qo+i-sim - s) ui^s, h{} 

s=l 

lo+l 



d 

jehe=+l s=l 



By functional independence of ee ^^^^ it follows that the function ob- 



e=±l 

lo+l 



tained as M7,io+i ~ (^io+i-s{\I\ — s) Ui ^ does not depend on any dynamical variable. 



s=l 



Hence it defines the coefficient Qq+i(|/|), thereby proving the r.h. to order Iq + 1. 
Finally, relation (O) immediately implies the result of the theorem. □ 



3.3 Comments 

An immediate consequence of this theorem is the existence of linear triangular rela- 
tions between the BCat Ruijsenaars- Schneider Hamiltonians and the classical Koornwinder- 
van Diejen Hamiltonians when w{qj) = v{2 qj) v{qj). 

The explicit coefficient have to be computed order by order since at this time no 
general recursion formula is available. As an example we have worked out the first 
two functions: 

Uk,i = Uk,i + c^{\K\) and Uk,2 = Kk,2 + c,{\K\ - 1)Uk,i + c^m) 

sinh7(2 \K\ + 1) 



with ci{\K\) 



sinh7 



1 /sinh7(2|fs:| - 1) sinh7(2|ir| + 1) sinh47|fs:| 
C2{\K\) = - 



2 \ sinh7 sinh7 sinh2 7 
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No obvious pattern appears yet. As a consequence, an algebraic interpretation of the 
Koornwinder-van Diejen Hamiltonians in connection with the canonical Hamiltonians 
is still lacking. 

More generally, the theorem implies that a hierarchy of Poisson-commuting Hamil- 
tonians with the generic form given is uniquely determinated by the giving of the fam- 
ily of f-functions and the first Hamiltonian, or equivalently the first "potential term" 
uk,i- In the Koornwinder-van Diejen case, this first Hamiltonian is given in [|l^ . 



We wish to end this section with a conjecture on the classical Koornwinder-van 
Diejen Hamiltonians with a general one-body potential chosen as in Jill (this time 
dropping the restriction to w{qj) = v{2 qj) v{qj)). They have not yet been constructed 
by a Lax formalism. We expect that the suitable Lax matrix for this hierarchy may 
be obtained by multiplying the BCat Lax matrix by a 2 + 1 diagonal matrix: 
L\Y = LDw, with D]y-- = W(gj) >V(— g^). This one-body potential function W has 
to be determinated by integrability conditions. In addition, we conjecture that, after 
some canonical transformation on the dynamical variables, the Hamiltonians 7i\^ = 
tr{L\y) will take the same form as the Koornwinder-van Diejen Hamiltonians (|T^ 
up to the change of Uk,p into some lAw^p- The theorem will then apply, thereby 
yielding the full connection between BCiv-type Ruijsenaars-Schneider potentials and 
the classical Koornwinder-van Diejen Hamiltonians. 
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